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Abstract 

We consider the effect of exact gluon kinematics in the virtual photon-gluon impact factor at small 
X. By comparing with fixed order DIS scheme splitting and coefficient functions, we show that the exact 
kinematics results match the fixed order results well at each order, which suggests that they allow for an 
accurate NLL analysis of proton structure functions. We also present, available for the first time, a;-space 
parameterisations of the NNLO DGLAP splitting functions in the DIS scheme, and also the longitudinal 
coefficients for neutral current scattering. 

The study of the proton structure functions at small x is of phenomenological importance, given the partonic 
centre of mass energies now accessible in collider experiments at HERA |2] and the forthcoming LHC. 
However, the coefficient and splitting functions relating the structure functions to the parton distributions 
contain logarithms in the Bjorken x variable. Although QCD fits at ncxt-to-lcading order (NLO) in as 
describe the data well, there is some evidence that a resummation of \ogl/x terms would improve the fits 
[3]. This is accomplished in principle via the BFKL equation an integral equation for the unintegrated 
gluon 4-point function /(x, k\, fc|) whose kernel is known at next-to-leading logarithmic (NLL) order [Sj. In 
deep inelastic scattering, the moments of the structure functions J-i{N,Q^) — Fi{x,Q^)dx are then 

given by the high energy factorisation formula 

TdQ',N)^asJ^ -^h{k^/Q^)f{N,e,Ql)gB{N,Ql), (1) 

where qb is the bare gluon distribution at momentum scale (5§, and the strong coupling as is fixed at LL 
order. The hi{k'^/Q'^) are the impact factors coupling the virtual photon to the gluon. At present these 
are known to LL order only [7| |S]. Thus a full NLL order small x analysis of the structure functions is 
not possible. In section 1 of this paper we show how to obtain the quark-gluon splitting and longitudinal 
gluon coefficient functions from the impact factors. In section 2, we discuss the DIS scheme coefficient and 
splitting functions at NNLO, needed to compare directly with the small x expansion, as these have not been 
presented before. In section 3 we compare the splitting and coefficient functions obtained from |^ with 
the complete NLO and NNLO results, showing that the inclusion of exact gluon kinematics in the photon 
impact factors gives a good approximation to higher order effects at small x, and quite possibly contains the 
dominant information in a NLL (or higher) order calculation. 



1 P^^^ C^g^ with exact gluon kinematics 

After solving the BFKL equation for f{N,k'^,QQ), the latter two factors in equation |^ combine to give 
a regularised unintegrated gluon distribution, and hi can naively be interpreted as the coefficient function 
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Figure 1: Diagrams contributing to the LL impact factor. Two further diagrams are obtained by reversing 
the direction of the quark loop. 



hnking the gluon density with the structure function. However, in the case of F2, the impact factor diverges 
as /Q^ 0. One can understand this given that at 0{a^), F2 is proportional to the quark singlet parton 
distribution with no gluon contribution. One does not expect to describe this nonperturbative dependence 
using perturbation theory, and thus the impact factor diverges. One must instead consider solving the 
evolution equation for F2, via the quantity: 



dT2iQ^,N) _ f°° dk 



ainQ2 -^h2ik'/Q')f{N,k\Qi)gB{N,Qi), (2) 

which serves to define the impact factor h2. In a general factorisation scheme, one loses the simple inter- 
pretation of h2 as the coefficient function relating the gluon distribution to the structure function. Instead 
it represents a mixture of the coefficient C2g and the anomalous dimension jqg. If one chooses to work in 
the DIS scheme P|, where F2 is given by the naive parton model expression to all orders, then /i2 can be 
interpreted directly as the quark gluon anomalous dimension "Y^g^^- From equation the longitudinal 
impact factor /il is identified with the coefficient function C^g^ and does not diverge due to the fact the the 
longitudinal structure function vanishes at O^a'g). It is convenient to perform a second Mellin transformation 
on the factorisation formulae to unravel the convolution in fc-space: 

Tl{i.N)^ dk\k^)-'-':FL(k^N) = hLh)G{j,N), (3) 
Jo 

and similarly for equation 1^)1. 

Diagrams contributing to the impact factor are shown in figure ^ One may introduce a Sudakov decompo- 
sition for the 4-momenta k,r: 

I ^ aq + l3p + l_i_] (4) 
k — — q' + Xgp + kj_, (5) 

where p is the proton 4-momentum (light-like if one ignores the proton mass), q' — q + xp a second light-like 
vector involving the Bjorken variable x — /{2p ■ q), and s = (p + q)^. The on-shell requirements for the 
intermediate quarks {{I + fc)^ — {q — 0^ — 0) then lead to the relation: 



(6) 



with A = a{l — a)A'^, — — g^, and = l± + {1 — a)k±. At LL order, the momentum fraction Xg of the 
incident proton carried by the gluon is undetermined, as equation ® implies the difference logXg — logx is 
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finite as a; — *■ 0. By imposing correct kinematics for the gluon, one includes in the impact factor significant 
higher order information. 

The resulting N dependent factors h2{'j,N) and hL{'-f,N) can be found in ^U], and from them one may 
derive estimates of 7^^"^ and C'f/'^ at fixed order in as- One first expands the relevant impact factor as 

a Taylor series in 7 with coefficients h^j"\ In solving the BFKL equation, 7 is identified as the anomalous 
dimension j{N) given at NLL accuracy in |S] (any further accuracy would require knowledge of the NNLL 
BFKL kernel). One thus has: 

00 

Cff('=)(a5,iV) = ^/.i")(7V)[7(as,iV)]("), (7) 

71=0 

as the exact kinematics result for the coefficient function up to NLL order, and similarly for ^^g^'^'^^ in terms 
of /i2- The BFKL anomalous dimension has a perturbative expansion: 

00 

7(as,iV) = ^aS/„(iV), (8) 

n=l 

so that equation Q is a power series in as, beginning at 0{as)- The corresponding j:-space expressions 
are given in Appendix A. One may compare order by order with the complete results for C^g and 7^3. The 
corresponding MS functions have been computed up to 0{ay) (111 1121 f?!! 1131 IT^ 1151 ITC) . However, for a 
direct comparison with the exact kinematics results one needs the corresponding results in the DIS scheme 
rather than the conventional MS scheme. 



2 DIS scheme splitting and coefficient functions 

The NNLO singlet and non-singlet splitting functions have recently been computed in the MS scheme |15II17| . 
along with the 0{a\) coefficient functions for neutral boson exchange ,16s TS|. They are extremely lengthy - 
for example, the typeset NNNLO gluon coefficient C2g is around twelve pages long ^S]- The NNLO splitting 
functions are shorter, but like the coefficients are made more complicated by the nature of the algebraic 
functions involved. All of the NNLO coefficient and splitting functions involve combinations of harmonic 
sums in A''-space, which after inverse Mellin transformation yield harmonic polylogarithms JHI in a;-space 
(up to weight five at this order). These are non-standard functions and thus must be generated numerically 
|20| . The combination of length and numerical complexity makes it is infeasible that the complete results can 
be immediately used in phenomenological applications. Instead one may parameterise the results in x-space 
in terms of simple algebraic functions, with a precision that far exceeds that due to higher order corrections. 
The Mellin transforms of the parameterisations then give suitably accurate A^-space representations. Pa- 
rameterisations of the MS functions are given in |15l 1171 1161 IT^ . From these we have derived corresponding 
representations of the DIS scheme quantities, accurate to within a percent apart from near the zeros. Our 
results are presented in Appendix B. 

In transforming between the MS and DIS schemes, we follow the argument presented in |21l I22j . The DIS 
scheme is characterised by the singlet structure function F2 having the same form as the naive parton model 
to all orders |5]: 

F2s[x, Q2) ^ tP'^ = C^®i:^ + C2g ® g^, (9) 

where E = X]i(9i + 9i) is the singlet quark density, and the factorisation and renormalisation scales have 
been chosen as Q^. The factorisation scheme independence of F2 then imposes a transformation between 
the DIS and MS scheme partons. There remains an ambiguity in the definition of the DIS gluon. However, 
the momentum sum rule fixes: ^ 

/ dxx[j:{x,Q^)+g{x,Q^)]^l (10) 
Jo 
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in both schemes. In MeUin space, this becomes: 

+gDiS(^) ^ S^^(iV) + .9^(Ar) (11) 
for = 1^. One may remove the ambiguity by extending equation pi|) to all N, and one obtains: 

To obtain the splitting functions, one differentiates equation H12|l with respect to and rearranges yielding: 

P""'' = ( Z + fiias)^) ® Z-\ (13) 



das J 

f P P \ 

where P = ( I . Substituting the perturbative expansions of the MS scheme coefficient and 

\ ^gq ^99 J 

splitting functions^, along with the QCD /3 function^, one can derive the DIS scheme results order by order 
in as- The explicit transformations at 0{a^g) are: 



p(2)DIS 

qq 


_ p(2)MS , 

qq ^ 




^ p(0)MS 




p(0)MS I 
99 ""^ 




> p(l)MS I 

^ gq ^ 


^ ^a)Ms ^ 


,p(l)MS 






p(0)MS^ 






p(0)MS^ 


^(1)MS 




5p(0)MS^ 








pirns ^ 








o ^(1)MS ^(1)MS „o ^(2)MS 






















p(2)DIS 

9q 


_ p(2)MS 

gq 




. p(0)MS 

qq 




p(0)MS 
99 




> p(0)MS , 
^ ff9 ^ 




?) p(0)MS 

^ gg 






p(l)MS 

qq 


^(l)MS 


^p(l)MS_ 




p(l)MS , 




J pa)Ms 










^(1)MS 
' ^2q 






^(1)MS 
•-"29 












p(l)MS^ 
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(14) 



(15) 



_ p(2)MS 
99 


^^2) MS ^ 


, p(0)MS 

^ qg 






p(0)MS g 
99 ^ 








p(0)MS g 
9ff ^ 
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gg - ^ gg ^2q ^ -^gg ^2g ^ -^gg "-"2g "A* qg ^2g "A" -^gg 



(16) 



p(2)DIS 

qg 


_ p(2)MS _ 
99 




p(0)MS 

^ qg 




p(0)MS 

gg 


cg^'''®n^°)'''+^2r'^ 


?) p(0)MS 
^ 93 






, p(l)MS , 

^ qg ^ 


^(l)MS 
~ L'2g 






p(l)MS , ^(1)MS p(l)MS 








~ p(0)MS g 

^ gq ^ 






pWMS^ 


^(l)MS_^(l)MS^p(0)MSg 








5p(0)MS^ 











(17) 



^This corresponds to our choice of Mellin variable and that of IIUI . The alternative definition f{N) = x^^ ^•'/(^) is also 
in common use, and in that case the second moment is constrained. 

^Conventionally, P^"^ is the coefficient of a"+^, where a = as/{47r); C'jj'^-j^ is the coefficient of a". 
^Here /3„ is the coefficient of a"+^. 
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Non-singlet quark combinations transform according to: 



DIS 



(18) 



which has the form of equation H12|) but with a trivial transformation matrix. Hence one obtains for the 
non-singlet splitting functions relevant to neutral and charged current scattering [Tzj : 



p+,-(2)DIS ^ p+,-(2)MS ^ p^d'^MS , 



'2q 



^(1)MS „ . ^+-(2)MS ^ ^(1)MS. 



'2q 



The pure singlet splitting function is given by: 

p(2)DIS ^ p(2)DIS _ p+(2)DIS 



(19) 
(20) 



The F2 coefficient functions are simply defined to all orders in the DIS scheme. For the longitudinal coeffi- 
cients, one considers: 



Pr — f r'MS r^MS r'MS 
Using the transformation equations p2|l and l|18(l . one finds: 

( cYf cYf GDIS )^(cYf cYf 



(21) 



Z 

a 



MS 
2ns 



(22) 



Explicit results at 0(ag) after substituting the expansions of the MS scheme coefficient functions are: 



^(3)DIS ^(3)MS 
^Lg — '-'Lg 














_ ^(2)MS ^ 
















^(l)MS 

•"is 

























(23) 



^(3)DIS 



^(3)MS 

^Lq - 






- 4?" ^ 






^(1)MS 


_ ^(2)MS ^ 
















(^(1)MS 
Lg 

























(24) 



^-K3)DIS _ ^(3)MS ^(2)MS ^(1)MS ^(1)MS ^(2)MS ^(1)MS ^(1)MS ^(1)MS 

•-^Lns — '-'Lns ~ Lns ^ '-'2ns ~ ^ Lns ^ '-'2ns + Lns ^ ^2ns ^ '-'2ns ■ K'^^) 

Then the pure singlet coefficient is given by: 

^(3)DIS _ ^(3)DIS ^+(3)DIS fr,f.^ 
^Lps — ^Lqq ~ ^Lns ' V^OJ 

The transformation terms were evaluated in iV-space, and divergent high and low N limits were then ex- 
tracted. Up to 0{1/N), one has a choice in how to extract the high N piece. We have chosen this in such a 
way as to lead to simple plus distributions and logarithms of (1 — x) in the x-space functions. The remain- 
ing finite functions as A'^ — > 0, oo were parameterised in cc-space by evaluating the inverse Mellin transform 
numerically. Finally the transformation terms were added to the existing MS parameterisations. Thus the 
plus distribution and small- a; divergent terms are exact up to truncation of the coefficients, as also are the 
parts of the log(l — x) terms not involving (1 — x) log(l — x). 
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The coefficients of (5(1 — x) in P,+ , Pgq and Pgg have been modified, and 6{1 — x) contributions added to Pqg 
that should in principle be absent. This, following refs. |15lll7[[THlll6| . is to increase the iV-space accuracy, 
such that the parameterised functions satisfy the momentum sum rules: 



(27) 
(28) 



for iV = 1 ^. One can also introduce such terms into the longitudinal coefficient functions, by fitting to 
numerical values of the low integer moments. We choose not to introduce these, however, given the size of 
these effects (no more than a few parts permille) do not exceed the uncertainty of the parameterisations. 
We have checked all of our expressions against known numerical moments |23| . 



Particularly noteworthy is the singularity structure of the DIS scheme functions as x 1 . One sees that the 
singlet quark splitting functions contain plus distributions up to 1^2 ^, or \og^{N) in Mellin space. However, 

(2)DIS 5 

Pgq contains more singular terms up to = log {N). One can understand this by considering what 
happens in the MS scheme. There \og{N) terms arise in the coefficients C2q and C2ns as a result of soft gluon 
emission from the quark probed by the virtual photon. As a; — > 1, there is insufficient phase space for the 
emission of real gluons, and thus an incomplete cancellation between singularities arising from virtual and 
real emission. The leading logarithms in N exponentiate |24[l25j . and the sub- leading logarithms can also be 
resummed [^Hl HZI 123 EH] • Combining the known resummation and fixed order results allows knowledge of 
the four leading towers of high N logarithms in C2q to all orders in as (SOI^- In the DIS scheme there are no 
such logarithms in the coefficients, as C2q is defined trivially to all orders. Instead the soft gluon resummation 
effects enter the splitting functions. The leading log(iV) terms in C|^| are produced by exponentiating those 
in 7^J®. This follows from the iV-space evolution equation for the non-singlet quark density: 



aiogQ2 



Ins W 



(29) 



which is easily solved to give: 

9°f(Q^) = gJ?f(Qg)exp 

= g°f (Q2)exp 



/ ns 



das 

W7) 



E\ ' / o (n)DIS m+n-1 , 

"S(y5) n=Om=0"'"s(Qo) 



(30) 



where the Cn,m are coefficients obtained after substituting in the perturbative expansions of the (3 function 
and anomalous dimension. Performing the integration in the exponent gives: 



(O)-DIS 



DIS 



DIS 



(Qo) 



/30 



exp ■ 



' (l)DIS 

Jns 



PI 



in 



in 



where the ellipsis denotes terms giving rise to sub-leading logarithms. Given that 7!"'^^^ 
Jns^^^ ^ ^og^{N) as iV — > 00, the leading logarithms in the exponent come from the term in 7, 
of the non-singlet structure function in the DIS and MS schemes is: 



(31) 

^ logiV and 
il\ The form 



DIS 

H'ns 



c. 



MS 
2ns 



(32) 



■^This also implies that the nj independent part of Pgg^^^'^ should vanish, given that Pqg has no term at 0{rJj). 



^See Appendix B for the definition of these functions. 
^This analysis has very recently been extended to include even higher order logarithmic corrections to DIS and Drell-Yan 
type processes 1311 . 
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Thus from equation (|30|l . ones sees that the leading powers of log(iV) in the MS scheme non-singlet coefficient 

function are generated by exponentiation of those in the DIS scheme NLO anomalous dimension jns^^^ (the 
LO anomalous dimension is independent of the factorisation scheme, and thus the prefactor in equation 
(|31|l is also found in the MS scheme). The next-to-leading log(iV) terms in are not so straightforward, 
but are determined by the exponentiation of a mixture of jnJ^^^ and jns'^^^ , and so on for the other 
sub-leading logarithms. A similar argument relates the leading \og{N) terms in with those in jqq^^^^ ■ 

This explains the absence of more singular logarithms ~ log^(iV) in the DIS scheme jq'^\ as the highest 
power of log(iV) is limited by the fact that it cannot exceed the power obtained by exponentiation of the 

leading log term in jqq^^^^ . Looking at equation (|14|l . the transformation terms in Pqf^^^ — > Pqq^^^^ involve 
the combination /3oC^J^^^ ® C^]^^^ - 2poC^f^^- Thus log'*(Ar) ter ms in C^^^ are cancelled by the com- 
bination [C'2j'''^^]^/2! due to the exponential structure of the leading logs in the coefficient function. 

The I?4 term in Pgq^^^^ corresponds to a next to leading high x divergence in 

C¥q^ ('^ allog^l^) in Mellin 

space), arising from the terms in equation ()15|l : 



p{2)DIS 

gq 



(2)MS 
2q 



p(0)MS 
99 



gg 



(1)MS 
2q 



c. 



(1)MS 
2q 



p(0)MS 
99 



•C. 



(1)MS 
2q 



(33) 



The small x limit of Pqg will be discussed in section 3 of this paper. Looking at the other splitting 
functions, one may verify the LL relations 0: 

^gq ~ Ca ^ 27r -RsJ ' > 

where Ca = 3, Cf = 4/3 are the QCD Casimir invariants and Tr — 1/2. These relations are also true at 
LL order in the MS scheme. 



The non-singlet and singlet splitting functions are plotted in figures [3 and 13 respectively. The singlet func- 
tions have been multiplied by x to alleviate the small x divergence. Aside from the differences at high x 
discussed above, one sees that the DIS scheme functions are more divergent at small x. This is analogous 
to the high x behaviour - in changing schemes one transfers divergences from the quark singlet coefficient to 
the splitting functions. Note also the qualitatively different structures at intermediate x in the two schemes. 
Each of the singlet splitting functions develops an extra turning point in the DIS scheme. 

The NNLO Pqg develops a negative dip in the DIS scheme at intermediate x, before increasing again as 
a; — > 0. Together with the large negative dip at high x this gives a negative result at intermediate x when 
convolved with a model gluon distribution which is more singular than the splitting function. We return 
to this feature in section 3. In fact, the qualitative structure of the DIS scheme splitting function can be 
reproduced from the truncated transformation: 

- Pqg^^ + ^2?^ ® P^g^^ - 2/3oCg)^. (35) 

The longitudinal quark and gluon coefficient functions are shown in figures ^ and [S] The two gluon 
coefficients are extremely similar. The pure singlet and gluon coefficients share the same small x limit, as 
the LL coefficients are the same in both schemes. There is, however, an extra turning point in the DIS 
scheme pure singlet function at higher x. Also of note is the negativity of the non-singlet coefficient at small 
X, a property also shared by the NLO result such that the complete non-singlet coefficient is negative at 
small X. However, it is not divergent as x — > so that convolution with a suitable non-singlet test function 
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Figure 2: The non-singlet splitting functions in the DIS (solid) and MS schemes (dashed), for nj = A. 



does not give a negative non-singlet structure function (see figure EJ. 

We now compare the DIS scheme quark-gluon anomalous dimension and longitudinal gluon coefficient with 
the corresponding results derived from exact gluon kinematics. 

3 Comparison of exact kinematics with NLO and NNLO results 

At 0{as), the exact kinematics results correspond exactly with the complete results as at this order, all the 
relevant diagrams are included in the impact factor calculation. The imposition of exact kinematics then 
supplements the x dependence that is missing when evaluating these diagrams in the LL limit. At higher 
orders, one can compare the complete TV-space functions with the estimates obtained from the modified 
impact factors. Expanding in N one finds: 



^WDisie) ^ Mmiy _ (102.97V + .2 1 40712.) + (172. 17^ + A39hj})N 



-{246.0nf + .659Snj)N'^ + O{N^); (36) 

^(2)ms(e) ^ f i_ L + (755571/ + 118.07i^ + .0 1 68271^) 

- (1508971/ + 204.071^ + .06 9 587i^)Ar + (25 1 6671/ + 3 1 2. Iri^ 

+ .U62n^f)N^ + 0{N^), (37) 
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Figure 4: The longitudinal quark coefficient functions in the DIS (solid) and MS schemes (dashed), for 
n/ = 4. The pure singlet coefficient has been multiplied by x. 
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Figure 5: The longitudinal gluon coefBcient function in the DIS (solid) and MS schemes (dashed) with 
n/ = 4, multiplied by x due to the divergence at small x. 
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Figure 6: C^nl^ *° ^i'^s) convolved with q„s = x °'^(1 — a:)^ (note an additional factor of a; is not included). 
Despite the negative sign of the coefficient at low x, the structure function FLns ~ xChns ® Qns is positive. 
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^ 20000 




Figure 7: The DIS scheme NLO and NNLO quark-gluon anomalous dimensions (sohd) (multiplied by x to 
remove the small x divergence) alongside the results obtained from the LL photon-gluon impact factor with 
exact gluon kinematics (dashed). Also shown are the LL approximations (dotted). 



The complete NLO and NNLO results give: 

5.072nj)N 

(38) 

399. 5n^ + .5926n^) 

(39) 

The leading logarithms in x (most divergent terms as — > 0) are correctly predicted from the resummation, 
and one sees that the next to leading terms in ^^g^^^ are well estimated by the exact kinematics expression 
(within 2% at NLO and 7% at NNLO, for n/ = 4). Accuracy falls off for higher order terms in A^, although 
these are not associated with small x divergence. The x-space functions are shown in figure [3 The exact 
kinematics results qualitatively reproduce the structures of the complete results, even at high x. They 
clearly do much better than the LL terms at approximating the splitting functions. Note that the small- a; 
behaviour does not set in until rather low x, as can be seen by the splitting function only turning positive 
for a; < 2 X 10~^ at NNLO (for nf ~ 4). The qualitative trend is that at higher order in as, the splitting 
function turns positive at lower x. We have confirmed, for example, that the NNNLO exact kinematics 
splitting function does not turn positive until x < 10"''. A good estimate for these values is obtained by 



34.67nj 
A^ 

+ (-374.6n/ + 11.70nj)N^ + 0{N^); 



7gg^^^^ = + (-109.371/ + .8889n^) + (233.6n/ 



^T'' - + -3165n, + 30.19n^, ^ ^^^^^^^^ 

+ (-34493n/ + 1589nJ - 6.121nJ)A^ 

+ (73141n/ - 4389n^ + 27.5Sn^f)N^ + 0{N^). 
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Figure 8: Example diagrams which contribute at 0{a^n'l) to the structure functions. 



approximating the exact kinematics sphtting functions by their asymptotic hmits as a; — > 0: 

a;P^'2)c/s(e) _ uiA7nf\og ^ - (2635n/ + 49.53n^) + . . . ; (40) 

^p{3)DiS{e) _ iiQYij^^i,^ iQg2 1 ^ (_78095n/ - 1410n^) log - + (248414n/ + 6924n^ + 8.2654) + . . . , 

(41) 

where the ellipses represent terms vanishing in this limit. Setting xPqg^'^^^'"'^^ — gives the approximate 
value a; = xo at which the LL terms begin to dominate over the sub-leading logarithms. For nj = A, one 
finds — 3 X 10"'^, 1 X 10^^ at NNLO, NNNLO respectively. The lower value of xq with increasing order 
of as implies that the leading small- a; resummation effects become less important phenomenologically at 
higher orders, as sub-leading logarithms dominate until very small x. 

From equations (|38fl and H39|l . we note that in the non- leading logarithmic terms, contributions involving 
higher powers of n/ are estimated poorly - including being of the wrong sign. This is expected given that 
higher powers of n/ in the perturbative contribution to the structure functions may arise from diagrams such 
as those shown in figure |S1 with fermion bubbles in the vertical rungs of the gluon ladder and in the quark 
loop at the top of the diagram. The former are included in the NLL BFKL anomalous dimension'', but the 
latter are missing in the exact kinematics calculation due the LL nature of the impact factor. However, one 
can see that the higher order n f terms do not constitute a very significant contribution relative to those at 
0{nf). Similar expansions for the longitudinal coefficient are: 

^(2)Dis(e) ^ _5^3337v ^ 22n; + .0329271^) + (62.52,1/ + ■U27n))N 

-(86.88n/-f .2551n^)7V2 + C'(7V3); (42) 

~(S)DIS(e^ 409.5nf -l2AQnf + 1.727 nl ^ „ „^ 

= + + (2127n/ + AQ.Un) + .01561n3) 

- (34 3671/ -f 68.9571^ -f m?>92n))N + (5345n/ -f 89.737i^ 

+ mZ2Qn))N'^ + 0{N^). (43) 

^Fermion bubbles in the bottom vertical rung of the ladder are not in the NLL anomalous dimension, but contribute to the 
scale of the coupling. See 1321 . 
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Figure 9: Diagram contributing a term cx flf^ to the structure function. 



The complete results give: 

C(^)DIS ^ „ 5.333n/ ^ ^_Q 229nf + .8889nJ) + (80.69n; - 4.850n^)7V 

+ (-133.8n/ + 10Mn})N^ + 0{N^); (44) 

~Ci\nr<^ 409 5nf -2076n <• + 102.4n? 

= + -j;^ + (4730n/ - 340.8n^ - .n39m,n} 

+ .5926n^) + (-9211n/ + 854.1n^ - .4340/;fin^ - 5.973n})N 

+ (20054n/ - 2251n^ + .0826Afll^nj + 25.7Anj)N'^ + 0{N^), (45) 

where flf^ —< e >^ / < >, taking averages over the active quark charges. The estimation of NLL terms 
is not as good as for jqg, even for the 0{nf) contribution. Again taking n/ = 4, the NLL term in the NNLO 
coefficient is estimated to within 35%. Nevertheless, the exact kinematics results are in good qualitative 
agreement with the complete results. The term in flfi will not be estimated by the exact kinematics 
calculation due to missing diagrams of the type shown in figure |3 Also, this term is not associated with a 
small X divergence at 0{a^g). Higher order terms in Uf and flfi are not very significant contributions. 

The x-space functions are shown in figure ITUI Again the exact kinematics results have the same qualitative 
behaviour as the complete results at both small and large x, whereas the LL approximations are compara- 
tively poor. 

The greater accuracy in •jgg can in part be attributed to the derivative of F2 in log(Q^). In Mellin space, 
this amounts to multiplication of the sum of the transverse and longitudinal impact factors by j{N), which 
suppresses the differences noted above by as- 

The a;-space functions will ultimately be convolved with parton distribution functions. Hence it is necessary 
to check the behaviour of the x-space exact kinematics expressions when convolved with a suitable gluon 
distribution. Following we convolve with the model gluon distributions: 

xg{x)^x-'^-^l^xf; (46) 
xg{x)^x'>-''{l-x)\ (47) 

where the former corresponds to a high scale (~ 30GeV^), and the latter reflects the fact that the gluon 
can be valence-like (or even negative at low x) at low ~ IGeV^ ,33. .341 . One expects resummation of 
small X terms to be more important at low Q^, due to the higher value of as. The results for P^g' and P^g' 
are shown in figures HI] and ini Results for CY;^"^' ° ® g &nd C^j^l"^' ® g are shown in figure 1131 and 1141 
For both the splitting and coefficient functions, the exact kinematics results qualitatively approximate the 
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Figure 10: The DIS scheme NLO and NNLO longitudinal coefficient functions (solid) (multiplied by x), 
compared with the results from the LL impact photon-gluon impact factor with exact gluon kinematics 
(dashed). The LL approximation (dotted) is also shown. 
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Figure 12: {Pqg^^^^ (x) 'S>g{x, Q^))/g{x, Q^) with nf = 4, showing NNLO (solid), exact kinematics (dashed) 
and LL (dotted) results. 
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Figure 13: Cl^ C 



)g with Uf = 4, showing NLO (sohd), exact kinematics (dashed) and LL (dotted) results. 



structure of the complete results. Comparing them with the LL terms at small x, after convolution with the 
gluon, one sees that they are much closer to the complete results. The exception is C'}^g , where the LL terms 
convolved with the gluon distribution are closer to the complete result at low x, aided by the fact that at this 
order the coefficient function has no next-to-leading small x divergence. At NNLO, where NLL terms are 
present, the exact kinematics results perform better at small x. The exact kinematics and complete results 
generally agree more at the lower momentum scale. This is due to the less singular gluon distribution at low 
Q^, and the small x part of the coefficient playing a more dominant role. However, at higher Q^, the effect 
of a more singular gluon distribution will be compensated in part by a lower value of as, and resummation 
becomes less important. 



The NNLO exact kinematics splitting function gives a negative result when convolved with the more singular 
gluon, turning positive only at very low x ~ 10~^. This can be attributed to the large negative dip in the 
exact kinematics function (see figure Cj) at intermediate x. Given that the gluon is more singular than the 
splitting function, the low x limit of the convolution is dominated by both high and low x information in 
the splitting function. To see how this works, consider the model splitting function: 



P= -+BS(l-x), 
x 



(48) 



where the first and second terms give the dominant behaviour at small and high x respectively. Consider 
convolving this with the following "gluon" : 



c/ = x°'e{xQ - x), 



(49) 



which is singular or valence-like at low x depending on whether a < or a > 0, and vanishes at high x. One 
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results. 
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then has: 



91og(Q2) 



dixf) 



xP® f 



A' 

B-- a;" + 



(50) 



a 



a 



If a > 0, the small x term in the splitting function dominates the convolution. If on the other hand a < 0, 
the bracketed term in equation H50|l gives the leading small x behaviour, which is a mixture of both the 
small and large x terms of the splitting function. This also accounts for the lack of a common small x limit 
in the left-hand plots of figures m El ^1 as each of the three splitting functions has a different high x 
behaviour. Note that the complete NNLO Pqg also has a negative dip at intermediate x. This leads to some 
negative behaviour after the convolution, but not for x < 0.05 in figure [T^ 

4 Conclusions 

We have shown that the imposition of exact gluon kinematics in the LL virtual photon-gluon impact factor 
gives a good approximation to the NLL parts of the NLO and NNLO splitting and coefficient functions 
Pqg'^^ and C^g^ up to 0(q;|). The qualitative behaviour is also good over the whole x range. We see this 
both by examining poles in iV-space and also convolving the x-space functions with suitable model gluon 
distributions. Hence in the absence of the full NLL impact factor ®, we have confidence that the exact 
kinematics results can be used for an accurate NLL analysis of the proton structure functions. 

It may also be possible to impose exact kinematics in the impact factors for heavy quark production |38| . 
In this case, however, one needs to define a suitable factorisation scheme in order to interpret the impact 
factors in terms of splitting and coefficient functions. 

The NNLO DIS scheme splitting and longitudinal coefficient functions (excluding the coefficients for charged 
current scattering) have been parameterised and presented here. There are significant qualitative differences 
with the MS scheme results, particularly in the appearance of divergent high x terms from soft gluon 
resummations in the splitting functions. These functions are available on request and can easily be applied 
for parton analyses in the DIS scheme at NNLO. 
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A Appendix: The exact kinematics splitting and coefficient func- 
tions 

The A'^-space splitting and coefficient functions derived from the exact kinematic impact factors involve the 
function ^^{N) = V{N) /V {N) and its derivatives, where r(A^) is the Euler gamma function. The -0 functions 
can be expressed as analytically continued harmonic sums |39| , which one can then inverse Mellin transform 
to a;-space. For brevity we define: 



Lo = log(x), Li = log(l - x). 



(51) 



Then the results are: 




(52) 



Calculation is, however, in progress 135113611371 . 
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^{3)Dis{e)^^^ = ([3468x - 21Q0x^]Ll + [-mx-^L\ + (2312a; - bUx'^ - ISOOa;^)^! - 2820 

- 11320x2 + 13432a; + (192Li2(a;) + 3277^ - 2176/3)a;"^]Lo + [SOOx^ + 256x-^ 

- 1156x]L? + [-1412 + 1800a; + (3277^ - 384 - 192Li2(a;))a;"i - 4a;"2]Li - 1032 

- 40892/9a;2 + 6952a; + [64Li2(.x) - 192Li3(l - a;) - 384Li3(a;) - 12388/9 
+384C(3) - 32/37r2]a;-i)n/ 

+ ([272/9a; - 80/3a;2]L2 + [(-160/9a;2 - 64/27a;-i + 544/27x)Li - 992/27a;2 
+ 64/9x - 16]Lo + [-272/27a; + 32/27a;-i + 80/9a;2]L2 + [i60/9a; - 64/9a;-^ 

- 304/27 + 16/27x-2]Li + 1168/27 + 5024/81x2 - 2800/27x + [-128/81 



a; 



+32/8l7r2-64/27Li2(x); 
+ 16/243 (a;(l - x)[iLl + 2Liio + 14^0 - i?] - Lq ~ [x^^ ^ i „ 2x\Li 
-2-7x^ + lQx~x-^)n); (53) 

p(i)Dis{e) ^ |^j92 i20x^ - U6x]Lo + 1048/3^^ + 44 + 104/3x-i - 384x) n/ 

+ ([16/27x(x - if + 8/27]Lo + 8/27 + 16/9x(x - 1)) n}; (54) 

Pg(2)'D^sW (x) ^ ([1350x2 - 1734x + 1587]io + [-mx-^Lj + (-1156x + 900x2 ^ -^Q^g _ 560x-i)Li 

- 890 + 10160x2 - 11340X + (192Li2(x) + 327r2 - 2272/3)x"i]Lo + [578x - 529 

- 450x2 + 280x-i]i2 + [706 - 900x + (327r2 - 288 - 192Li2(x))x"i - 2x"2]2,^ 

+ 5992 + 86146/9x2 - 12248x + [136/37r2 + 384C(3) - 272Li2(x) - 29842/9 - 384Li3(x) 
-192Li3(l - x)]x"i) Uf 

+ ([40/3x2 + 92/9 - 136/9x]ig + [(184/27 - 64/27x"i - 272/27x + 80/9x2)Li 

+ 1696/27x2 - 1448/27 - 1184/27x]io + [-92/27 + 32/27x^1 + 136/27x - 40/9x2]l2 

+ [8/27x^2 _,_ 152/27 - 32/9x-i - 80/9x]Li + 800/27 - 472/81x2 + 88/3x 

+ [32/8l7r2 - 4304/81 - 64/27Li2(x)]x~i) nj 

+ 4/729 (2x(x - l)[3Ll + 2LoLi - Lf] + 3Lg + 2Li + 2[3 - 22x + 24x2]Lo 

-LI + 2[1 - 2x - x"2]Li + 50x2 _ 5g2. + g _ 2x"^) n^f, (55) 

where Li„(x) is the n**^ polylogarithm function. 

B Appendix: The DIS scheme spHtting and coefficient functions 

Here we present parameterisations of the coefficient and sphtting functions at NNLO in the DIS scheme 
12?. For completeness, all singlet and non-singlet splitting functions are given. The longitudinal coefficient 
functions are given only for neutral current structure functions, as the MS scheme coefficient functions for 
charged current scattering have yet to be published. First we define: 

A.= ~ , io-log(x), Li = log(l-x) (56) 

1 — X 
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Then the results are: 

P+mois ^ 785.o6Do - 2974.4Di + 645.33I>2 + 14669.3758^(1 - x) + 1868.3 + 6601.3a; + 243.6a;^ 

- 522.1a;^ + 77.391i? + [-2771.56a; + 3059.61]i^ + [2695.85a; + 13750]ii 

+ [1.5802 - 15.818.r]L(5 + 83.639L^ + [83.48ii + 915.18]Lg + [-272.00a; + 2497.50 

- 750.9^1 + 8314.3L2 + 544.00/(1 - x)]Lo 

+ Uf (-325.I8D0 + 403.89I?! - 78.222X>2 - 2150.5868(5(1 - x) + 12.951 + 217.65a; 
+ 358.28a;2 + 44.79a:3 + .95867.tL;5 + [2.573a; - 5M36]LI + [-118.68 + 10.503Li]Lg 
+ [-67.556/(1 - a;) - 155.46Li - 503.89L? + 33.778a; - 327.76]Lo - 4.6904L? 
+ [-152.43 + 167.97x]Ll + [-180.68a; - 1417.78]Li) 

+ rif (7.6750Do - 11.457Di + 2.3704D2 + 63.63585(1 - a;) - 4.8837 - 28.501a; 

- 17.293a;^ - .24667a;L^ + [1.1852a;/(l - x) - .59259a; + 3.5556]Lo + [11-457 
+3.9506a;/(l - x) - 4.3457a; + 10.817ii]io - 2.0000i? + 33.160ii) ; (57) 

PnP'^'^ - 785.06D0 - 2974.4Di + 645.33X»2 + 146595(1 - a;) - 42.670 + 10704a; + 297.0a;^ 

- 433.2a;^ + 1.4321L^ + 106.84L^ + [994.40 - 860.64ii]i^ + [-272.00a; - 630.82ii 
+ 2107.4 + 9310.9L? + 544.00/(1 - x)]Lo + 65.291L? + [3085.1 - 2771.5a;]i? 

+ [14503. + 2695.9a;]Li 

+ Uf (-325.I8P0 + 403.89151 - 78.222X>2 - 2150.0(5(1 - a;) + 75.786 + 413.96a; + 77.89.r2 
+ 34.76a;3 - [1.136a; + 7.4805]L^ + [.59212Li - 125.14]L^ + [-381.08Li - 67.556/(1 - a;) 
+ 33.778a; - 564.29L? - 321.26]io - 3.9570i? + [167.97a; - 149.42]L? 
-[1421.1 + 180.68a;]ii) 

+ n} (7.6750Do - 11.457Di + 2.3704D2 + 63.5855(1 - x) - 2.0572 - 55.288a; + [3.5846 

- .59259a; + 4.0479Li + 1.1852a;/(l - x)]Ll + [23.959Li + 12.039 + 3.9506a;/(l - x) 
-4.3457a;]Lo - 2.2760L? + 30.600Li) ; (58) 

P^f°'^ ^Hf (-193299 - 672088a; + 104121a;2 + 964027.t-'^ - 201675a;^ - 1327.613;"^ - 820.836L^ 

+ [-61102. 6L1 - 9741.91]l2 + [307888Li - 190993L? - 75017.3 - 196.207a;-i]Lo + 2332.86Xi 
-1876643L? + 10385.8iJ - 1121.25a;Li + 1876481a;i? - 10509. 5.TiJ + 24.88888^^) 
+ nj (530.035 + 7303.79a; - 937.737a;2 - 7925.90a;3 + 1105.83a;'' + 9.365933;"^ - 1.65094L^ 
+ 1.51450Li^ + [251.096Li - 31.8095]Lo + [-1502.08ii + 123.563 + 13313.6Li]Lo + 2649.16Li 
+13574.4i? - 2639.02a;ii - 13583.3a;i? + 1.18519if) ; (59) 
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pp)r>js ^ 2643.5D0 + 4425.87395(1 - x) - 20852 + 3968a; - 3363a;^ + 4848a;^ + 14214a;-^ - 144L^ 
+ 72Ll + [8757Li - 7471]L^ + [274.4 + 2675.8a;-i + 7305Li]Lo + 3589ii 
+ nf (-412.172X>o - 534.1666(5(1 - x) + 94680.9 + 423522a; - 62541. Olx^ - 569436x^ 
+ 120946a;'' + 1149.99a;-i + 18.9631L;^ + 660.814i;5 + [24297.9ii + 5133.55]L§ 
+ [1099250L? - 175012Li + 220.737x-i + 40461.3]Lo - 24.8889L^ + [2062.11a; - 1913.21]LJ 
+ [- 1093524a; + 1093454]i? + [-22404.4a; + 22442.5]Li) 

+ n} (-1.77778D0 + 6.44153(5(1 - x) - 19903.1 - 81663.3a; + 11472.20;^ + 114322a;^ - 24596.2a;* 
- 17.8171a;-^ - 81.0657Lg + [-5570.25Li - 1006.77]Lo + [-7493.22 - 215788L? + (85.25a; 
+37019.4)Li]Lo + [784.390 - 787.057a;]if + [212770a; - 212747]Lf + [162.579 - 283.399a;]Li) 
+ (14.399 + 15.108a; - 104.84^^ + 41.797a;^ + 33.545a;* + .44376L^ + [-22.307ii 
+ 2.6393]Lg + [-139.31Lf - 112.57ii + 9.5276]Lo - .26473(a; - l)i? + 140.68(a; - l)Lf 
+112.93(a;- l)ii); (60) 

P^Dis _ (5.2833(5(1 - x) - 396354 - 1679228a; + 400583a;2 + 2086958a;3 - 413461a;* 
- 3164.7.T"i + 19.852L;^ + [-1610.1 - 252.5a;]Li^ + [-130777Li - 18993]L^ 
+ [-152378 + 612746Li - 441.47a;-i - 3977765L?]Lo + 36.004^^ + [2923.9 - 3143.9a;]L? 
+ [4014417a; - 4013999]L? + [325972x - 323111]Li) 

+ n} (-0.2404(5(1 - a;) - 2697.2 - 16699a; + 11510a;^ + 8064.9a;^ + 30.195a;-^ - 10.237L^ 
+ [-32.113 + 11.70a;]i^ + [-1088.9Li - 427.67 - 98.07a;]L^ + [2782.9ii - 1477.3 
-28176L?]Lo + 2.6670L? + [29629a; - 29652]i^ + [7440.6a; - 7360.6]ii) 

+ (0.00135(1 -x) + 143.623 - 835.290a; + 540.973a;^ + 150.487a;^ + .326927i^ + 3.75075L^ 
+ [21.1805Li + 25.3458]i^ + [95.2390 + 196.232ii - 865.734L?]Lo + 881.848(a; - 1)LI 
+32.8575(a; - l)Li) ; (61) 

P^Dis _ -10172. 599Po + 2619.956X'i + 3026.479I?2 - 75.85220^3 - 118.5187^4 

- 17666.5673(5(1 - x) - 63856.3 - 226976;r + 5645.05.t2 + 370971a;-'^ - 80184.6a;* + 6133.90a;-i 

- 52.9383L* - 269.675Lg + [-972.9a; - 16883.4Li - 6887.42]L2 + [-25459.0 + 1189.3a;"i 

- 692833Lf + 122913Li]Lo + 89.4494L^ + [7443.02 - 7918.33a;]L? + [-666832 + 658887a;]L? 
+ [16626.4 - 32060.8a;]Li 

+ n/ (935.1848D0 - 550.4791Di - 56.29637D2 + 25.28401I>3 + 2589.95315(1 - x) + 35445 + 73884a; 

- 11203a;2 - 127979^^ + 27317a;* + 350.55a;-^ + 4.7407i^ + 312.26ig + [9521. lii + 108.6a; 

+ 2357.5]L^ + [99.282a;-i - 45381Li + 16599 + 260063L?]Lo - 14.2221.^ + [-1762.4 
+1847.8a;]L? + [-254267.T + 254637]i? + [-4256.9 + 4324.8a;]Li) 

+ nj (-12.698P0 + 17.185Pi - 3.5556P2 - 93.6748(5(1 - a;) - 103.10 + 809.52a; - 655.80a;2 

- 5.0491a;-i - SmUOLf, + [31.758Li - 20.430]i^ + [-273.77ii + 317.30L? - 144.58]Lo 
-3.5556L? + [-230.63a; + 241.15 + 3.5556a;"^]L? + [96.502a; + 11.852a;"^ - 110.71]Li) ; (62) 
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CiHt^^^i^) - -3634.5 + 5025.2a; - 614.77a;^ - 996.21x'' + (1 - a;)[8452.3Li + 4090.2L? + 175. 59L? 

+ 225.30Lf] - 3280.3LoLi - 1082.7LgLi - 911.45Lo - 81.823L^ - .72047^;] - 1780.0^0^? 
- 4059.2Li + 125.02LJ + 21.113L? + 1.6059Lf 

+ n/ (617.05 - 1670.5X + (1 - a:)[23.584Li - 106.82L2] + 1717.8Lo-^i 
+465.96LgLi + 171.90Lo + 6.9942Lg + 370.25Li - 45.190L?) 

+ nj (-17.038 + 35.968a; + (1 - a;)[-22.215Li + 23.829L?] - 66.179Loii - 12.884L^Li 
-5.I888L0 - .025315Lg + 24.794Loi? - 15.012Li + 2.3704L?) 
+ /^n"/ ([107.0 + 321. 05x - 54.62a;2](l - x) - 26.717 + 9.773Lo 

+ [363.8 + 68.32Lo]a;Lo - 320/81Lg[2 + Lq]) x; (63) 

Cf^f^^ix) ~ Uf (1769.7 - 441. 62x - [182.00Lo + 899.64]a;-i - (1 - a;)[23.584Li - 106. 82L?] 
+ 53648LoLi + 116041,2^^ _ 894.81Lo + 105.36Lg + (1 - a;)[81652Li + 3880. 7Lj] 
-76.310Lij - 8700.9LoL?) 

+ nj (4087.2 - 4143.1a; + 47.29x-i + (1 - a;)[2293.9Li + 654.38Lf] + 978.21Loii 
+ 2199.0LgLi + 1484.3Lo + 176. 52L^ + 18.3272.;] + 511.80Loi? 
-(1 - a;)[-22.215Li + 23.829LJ]) 

+ f^ilnf ([107.0 + 321.05a; - 54.62a;2](l - a;) - 26.717 + 9.773Lo 

+ [363.8 + 68.32Lo]a;Lo - 320/81L^[2 + Lq]) a;; (64) 

Cf^°^^{x) ~ n/ (-4573.1 + 77228a; - 70637a;3 - [409.506Lo + 2076.4]a;~i 

+ (1 - a;)[-8666.9Li + 267612L? - 4500. IL?] - 8146.1LoLi + 4257.5LgLi 

-4277.2Lo - 241.08L^ - 246.51^1] + 272818Loi? + .32800Li) 

+ n} (8878.1 - 14399a; + 5430.1a;^ + 102.40a;"^ + (1 - a;)[5143.2Li - 83.489Li] 

+7051.6Loii + 5593.0L^Li + 3258. 9Lo + 481.19Lg + 68.03411 + 516.40Loi?) 

+ n} (-287.66 - 494.46x + 782. 72a;^ + (1 - a;)[-614.31Li - 1547.3L?] - 32.680Loii 

-112.24L2Li - 132.42Lo - 26.899L^ - 2.6004Lg - 1490.4LoL2) 

+ flf^ji} ([-0.0105L? + 1.550L? + 19.72xLi - 66.745a; + 0.615a;2](l - a;) 

+ 20/27a;L^ + [280/81 + 2.260x]xLI - [15.40 - 2.201a;]a;L^ 

-[71.66 -0.121a;]a;Lo)- (65) 
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